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Equivalent-Continuum Modeling of Nano-Structured Materials 

Gregory M. Odegard . Thomas S. Gates, Lee M. Nicholson \ and Kristopher E. Wise 
NASA Langley Research Center, Hampton, Virginia 


Abstract 

A method has been developed for modeling structure-property relationships of nano-structured 
materials. This method serves as a link between computational chemistry and solid mechanics 
by substituting discrete molecular structures with an equivalent-continuum model. It has been 
shown that this substitution may be accomplished by equating the vibrational potential energy of 
a nano-structured material with the strain energy of representative truss and continuum models. 
As an important example with direct application to the development and characterization of 
single-walled carbon nanotubes, the model has been applied to determine the effective 
continuum geometry of a graphene sheet. A representative volume element of the equivalent- 
continuum model has been developed with an effective thickness. This effective thickness has 
been shown to be similar to, but slightly smaller than, the interatomic spacing of graphite. 

KEY WORDS: nanotechnology, nanotubes, continuum mechanics, molecular mechanics, finite 
elements 

Nomenclature 

Graphite orientation 
ap - Unit-cell base axes 

D - Graphite orientation indicator 

«p - Orientation integers 

Force field 

ET‘ - Non-bonded electrostatic potential energy 
E? - Vibrational potential energy of graphene sheet 
E" - Potential energy of nano-structured material 
E' dl 1 - Non-bonded van der Waals potential energy 

£° - Bond-angle variation potential energy 

E? - Bond stretching potential energy 
E l - Bond torsion potential energy 

EN - Bond inversion potential energy 

K ( p - Constant associated with C-C bond stretching in graphite 

K ; e - Constant associated with C-C-C bond-angle variation in graphite 

0, - Deformed bond-angle of C-C-C bonds in graphite 
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0 ; - Equilibrium bond-angle of C-C-C bonds in graphite 

p, - Deformed bond length of C-C bonds in graphite 
P ( - Equilibrium bond length of C-C bonds in graphite 

Truss and continuum models 

a - Elastic rod type of outer portion of truss representative volume element 
A' - Cross-sectional area of rod i of truss member type j 
A ni - Cross sectional area of carbon nanotube 

b - Elastic rod type of inner portion of truss representative volume element 
rf - Deformed distance between joints of rod i of truss member type j 

R ' - Undeformed distance between joints of rod i of truss member type / 

/ - Wall thickness of continuum nanotube 

E 8 - Young’s modulus of graphene sheet 

Y/ - Young's modulus of rod / of truss member type./ 

E" 1 - Young's modulus of carbon nanotube 

r - Mid-plane radius of continuum nanotube wall 

r” - Inner radius of continuum nanotube wall 

r m - Outer radius of continuum nanotoube w’all 

A c - Mechanical strain energy of the continuum model 

A 1 - Mechanical strain energy of the truss model 

v g - Poison’s ratio of graphene sheet 

Boundary conditions 

B\..B\{) - Displacement field constants 

- Vector of periocity of representative volume element face a 

f k - RVE and macroscopic coordinate system offset vector 
h - Height of macroscopic graphene plate 

u ( k a) - Displacement vector of representative volume element face a 

u k - Global displacement vector 

vi’ - Width of macroscopic graphene plate 

x k - Cartesian coordinate system of the representative volume element 

x k - Global coordinate system 

a - Representative volume element face number 
A - Prescribed displacement applied to macroscopic plate 

£w - Strain tensor 

Ou - Stress tensor 
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1. Introduction 


Nano-structured materials have excited considerable interest in the materials research community 
in the last few years partly due to their potentially remarkable mechanical properties [1]. In 
particular, materials such as carbon nanotubes, nanotube and nanoparticle-reinforced polymers 
and metals, and nano-layered materials have shown considerable promise. For example, carbon 
nanotubes could potentially have a Young's modulus as high as 1 TPa and a tensile strength 
approaching 100 GPa. The design and fabrication of these materials are performed on the 
nanometer scale with the ultimate goal to obtain highly desirable macroscopic properties. 

One of the fundamental issues that needs to be addressed in modeling macroscopic mechanical 
behavior of nano-structured materials based on molecular structure is the large difference in time 
and length scales. On the opposite ends of the time and length scale spectrum are computational 
chemistry and solid mechanics, each of which consists of highly developed and reliable 
modeling methods. Computational chemistry’ models predict molecular properties based on 
known quantum interactions, and computational solid mechanics models predict the macroscopic 
mechanical behavior of materials idealized as continuous media based on known bulk material 
properties. However, a corresponding model does not exist in the intermediate time and length 
scale range. If a hierarchical approach is used to model the macroscopic behavior of nano- 
structured materials, then a methodology must be developed to link the molecular structure and 
macroscopic properties. Even though there is a long history of modeling bulk properties of 
materials based on molecular properties, a simple link between the firmly established disciplines 
of computational chemistry and solid mechanics has not been established. 

In this paper, a methodology for linking computational chemistry and solid mechanics models 
has been developed. This tool allows molecular properties of nano-structured materials obtained 
through molecular mechanics models to be directly used in determining the corresponding bulk 
properties of the material at the macroscopic scale. The advantages of the proposed method are 
its simplicity and direct connection with computational chemistry and solid mechanics. 

This approach consists of two major steps. First, the chemical structure and vibrational potential 
energy function of the nano-structured material, which is determined from computational 
chemistry, is used to obtain an equivalent mechanical pin-jointed truss model. Second, the truss 
representation of the nano-structured material is substituted with an equivalent-continuum 
model. The mechanical behavior of the continuum model closely approximates that of the nano- 
structured material. 

As an important example, with direct application to the development and characterization of 
single-walled carbon nanotubes (SWNT), the proposed method has been used to determine the 
effective geometry of a graphene sheet. A representative volume element (RVE) of the graphene 
layer has been modeled as a continuous plate with an effective thickness that has been 
determined from the nano-structured properties of graphite by using the proposed method. 
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2. Background and basic concepts 

A discussion of some of the necessary background and concepts associated with characterization 
and modeling of carbon nanotubes is required for the development of the proposed method. 

2.1 Carbon nanotubes 

Carbon nanotubes have become a primary focus in nanotechnology research due to their 
exceptionally high stiffness and strength. One of the fundamental issues scientists are 
confronting is the characterization of the mechanical behavior and properties of individual 
carbon nanotubes. Many experimental [2-6] and theoretical [7-11] studies have been performed 
on single- and multi-walled carbon nanotubes. In particular, deformation modes and overall tube 
stiffnesses have been closely examined. 

In order experimentally characterize the mechanical behavior of nanotubes, a thorough 
understanding of the physical properties, such as effective cross-sectional area and moment of 
inertia, and materials properties, such as Young's modulus and Poisson’s ratio, is necessary'. 
These quantities are traditionally associated with the macroscopic scale, where the characteristic 
dimensions of a continuum solid are well defined. The determination of these properties has 
been attempted in many of the studies cited above without proper regard to an acceptable 
definition of the nanotube geometry. Accurate values of macroscopic physical and mechanical 
properties are crucial in establishing a meaningful link between nanotube properties and the 
properties of larger structures, such as nanotube-reinforced polymer composites. Therefore, 
caution should be used when applying continuum-type properties to nano-structured materials. 

In many studies, it has been assumed that the nanotube “wall thickness” is merely the interatomic 
spacing of two or more graphene sheets [3, 6-11], which is about 0.34 nm in single-crystal 
graphite. While this simple idealization appears to have intuitive merit, it does not necessarily 
reflect the effective thickness that is representative of continuum properties. In order to avoid 
this problem, Hernandez et al. [8] proposed the use of a specific Young’s modulus, i.e„ Young's 
modulus per unit thickness. Even though this approach is convenient for studies concerned with 
the relative stiffnesses of nanotubes, it is of little use when modeling a nanotube as a continuum 
structure. Another proposed solution to this dilemma is to assume that the nanotube is a solid 
cylinder [12-14]. This method is certainly convenient, however, significant inconsistencies arise 
when comparing moduli data of single wall nanotubes (SWNT) and multi-walled nanotubes 
(MWNT) when both are assumed to be a solid cylinder. 

A continuum modeling approach, which relies on an effective wall thickness, was proposed by 
Yakobson et al. [15]. They developed an elastic continuous shell model to represent a carbon 
nanotube. In their model, a wall thickness (0.066 nm) was calculated based on an assumed value 
of Young's modulus of the nanotube (5.5 TPa). This approach is extremely useful in 
deformation analysis of nanotubes (e.g., buckling), however, it is not convenient when 
attempting to calculate the modulus of the nanotube based on empirical or theoretical data where 
the thickness is not known a priori. Govindjee and Sackman [14] utilized the Bernoulli-Euler 
beam theory to show' that for a MWNT in bending, a large number of atomic layers are necessary 
in order to be able to assume that the cross-section of the nanotube is a continuum. This type of 
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approach assumes that the nanotube is a solid cylinder instead of using a more realistic model of 
single or multiple hollow cylinders with inter-wall interactions. Ru [16] used an elastic shell 
continuum model to investigate buckling of double-walled nanotubes subjected to axial 
compression. It was shown that the van der Waals forces between the layers did not increase the 
critical axial strains that correspond to buckling. Specific thicknesses for the inner and outer 
walls were not discussed in detail. Ru [17] proposed that the effective bending stiffness of 
SWNT should be regarded as an independent material parameter not related to the bending 
stiffness associated with the Bemoulli-Euler beam bending theory. This simplification was 
suggested in order to explain inconsistencies associated with the use of the classical bending 
formula for SWNT using a wall thickness of 0.34 nm. 

It has been postulated by Cornwell and Wille [12] and Robertson el al. [1 8] that the strain energy 
induced by bending a flat graphene sheet into a circular nanotube causes an increase in the 
apparent stiffness of the tube as a function of the nanotube radius. Conversely. Halicioglu [19] 
has shown that significant radial and circumferential stresses are present in the nanotube that are 
also dependent on the tube radius. Therefore, it follows that in order to properly model the 
mechanical behavior of a SWNT using continuum mechanics, the effective geometry must be 
known. If the nanotube is modeled as a continuous hollow cylinder with an effective wall 
thickness, then a simple first step is to model the flat graphene sheet in order to determine the 
effective wall thickness. In the current study, the effective thickness will be calculated as an 
example of the proposed approach. 

2.2 Notation for graphite orientation 

In order to clarify some of the notation used in this paper, a brief description of the commonly 
used notation for graphite is necessary. A graphene sheet is composed of bonded carbon atoms 
arranged in a repeating array of regular hexagons (Figure 1). The simplest way of specifying the 
orientation of a graphite structure is in terms of a vector D that joins two equivalent points 
(atoms) of the graphene sheet. The vector may be expressed as: 

D = /7,j a p ( 1 ) 

where the index (3 = 1,2; ap are the unit-cell base vectors of the graphene sheet (Figure 1); «p is a 
set of integers where and the summation convention is used for repeated indices. It is 

common to express equation ( 1 ) in terms of the integers as the ordered pair {n\,m). 

2.3 Representative volume element 

To reduce the time and number of computations associated with modeling the graphene sheet, a 
representative volume element (RVE) for graphene was used in this study (Figure 2). The 
selected RVE allows each degree of freedom of the carbon atom associated with bond stretching 
and bond-angle variation in the hexagonal ring to be completely modeled by truss and continuum 
finite element model nodal-displacement degrees of freedom. Also, this RVE allows the 
displacements on the boundary of the proposed chemical, truss, and continuum models to 
correspond exactly. Furthermore, macroscopic loading conditions applied to a continuous 
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graphene plate can be easily reduced to periodic boundary conditions that are applied to the 
RVE. 

3. Modeling procedure 

The proposed method of modeling nano-structured materials with an equivalent-continuum is 
outlined below. The approach uses the energy terms that are used in molecular mechanics 
modeling for the development of a continuum solid. Therefore, a brief description of molecular 
mechanics is given first followed by an outline of the equivalent-truss and equivalent-continuum 
model development. 

3.1 Molecular mechanics 

An important component in molecular mechanics calculations of the nano-structure of a material 
is the description of the forces between individual atoms. This description is characterized by a 
force field. In the most general form, the total potential energy of the force field for a nano- 
structured material is described by the sum of many individual energy contributions: 

E" = E P + E* + E X + E W + E' ,m + E cl (2) 

where E?, E 6 . £ T , and £"' are the energies associated with bond stretching, angle variation, 
torsion, and inversion, respectively. The non-bonded interaction energies consist of van der 
Waals. E vdW , and electrostatic, E’ 1 , terms. Various functional forms may be used for these energy 
terms depending on the particular material and loading conditions considered [20]. Obtaining 
accurate parameters for a force field amounts to fitting a set of experimental or calculated data to 
the assumed functional form, specifically, the force constants and equilibrium structure. In 
situations where experimental data are either unavailable or very difficult to measure, quantum 
mechanical calculations can be a critical source of information for defining the force field. 

3.2 Truss model 

Due to the nature of the molecular force field, a pin-jointed truss model may be used to represent 
the energies given by equation (2) where each truss member represents the forces between two 
atoms. Therefore, a truss model allows the mechanical behavior of the nano-structured system to 
be accurately modeled in terms of displacements of the atoms. This mechanical representation 
of the lattice behavior then serves as an intermediate step between linking the vibrational 
potential with an equivalent-continuum model. In the truss model, each truss element 
corresponds to a chemical bond or a significant non-bonded interaction. The stretching potential 
of each bond corresponds with the stretching of the corresponding truss element. Traditionally, 
atoms in a lattice have been viewed as masses that are held in place with atomic forces that 
resemble elastic springs [21]. Therefore, bending of truss elements is not needed to simulate the 
chemical bonds, and it is assumed that each truss joint is pinned, not fixed. 
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The mechanical strain energy, A 1 , of the truss model is expressed in the form: 


' XX (>;' * i 


(3) 


where A/ and Y/ are the cross-sectional area and Young's modulus of rod / of truss member type 
j. respectively. The term r' - R' is the stretching of rod / of truss member type./, where R' and 
r 1 are the undeformed and deformed lengths of the truss elements, respectively. For equation 
(3), as well as all equations up through section 4.2, the summation convention for repeated 
indices is not used. 

In order to represent the chemical behavior with the truss model, equation (3) must be equated 
with equation (2) in a physically meaningful manner. Each of the two equations are sums of 
energies for particular degrees of freedom. The main difficulty in the substitution is specifying 
equation (3). which has stretching terms only, for equation (2), which also has bond-angle 
variance and torsion terms. No generalization can be made for overcoming this difficulty for 
every nano-structured system. A feasible solution must be determined for a specific nano- 
structured material depending on the geometry, loading conditions, and degree of accuracy 
sought in the model. 

3.3 Equivalent-continuum model 

For many years, researchers have developed methods of modeling large-area truss structures 
with equivalent-continuum models [22-27], These studies indicate that various methods and 
assumptions have been employed in which equivalent-continuum models have been developed 
that adequately represent truss structures. In this study, the truss and continuum models are 
assumed to be equivalent under the following conditions: 

1. Truss lattices with pinned joints can be modeled accurately with an equivalent-continuum 
model that is based on classical continuum mechanics. For this case, micropolar [28] 
continuum assumptions are not necessary. 

2. The models have the same degrees of freedom. 

3. The displacements along the edges of the RVE are identical for the two models subjected 
to the same static loading conditions. 

4. The same amount of thermoelastic strain energy is stored in the two models when 
deformed by identical static loading conditions. 

If these criteria are satisfied, then a particular continuum model, such as a beam, plate, shell, or 
three-dimensional solid, may be directly substituted for a discrete truss lattice. The parameters 
of the solid, such as the elastic properties and geometry, are determined based on the above 
criteria. In some cases the strain energy of the continuum, A c . can be easily formulated 
analytically and compared directly with equation (3) to obtain the equivalent-continuum 
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properties. In other cases, especially with complex geometries and deformations, numerical 
tools need to be used to determine the continuum parameters. 

Once an equivalent-continuum model has been determined, the mechanical behavior of larger 
structures made of the nano-structured material may be predicted using the standard methods of 
continuum mechanics. 

4. Example: effective geometry of a graphene sheet 

In this section, a graphene sheet is modeled as a continuous plate with a finite thickness that 
represents the effective thickness for the determination of continuum-type mechanical and 
physical properties. By using the methodology described above, the molecular mechanics model 
is substituted with a truss model and subsequently an equivalent-plate model. The continuum 
model may then be used in further solid mechanics-based analyses of SWNT. 

4.1 Molecular mechanics model 

The force constants used in this example were taken from the MM3 force field of Allinger and 
coworkers [29-31]. Due to the nature of the material and loading conditions in the present study, 
only the bond stretching and bond-angle variation parameters were used. Torsion, inversion, and 
non-bonded interactions were assumed to be negligible for the case of a graphene lattice 
subjected to small deformations. For this example, the vibrational potential energy of a graphene 
sheet with carbon-to-carbon bonds is expressed as a sum of simple harmonic functions: 

£“ =£*'(ft- p .f +IX(e,~e,) J (4) 


where the terms P, and©, refer to the undeformed interatomic distance of bond i and the 
undeformed bond-angle /, respectively. The quantities pj and 0i are the distance and bond-angle 
after stretching and angle variance, respectively (see Figure 3). K , p and K® are the force 

constants associated with the stretching and angle variance, respectively, of the chemical bonds. 
Using the parameters for the MM3 force field [29-31], the force constants used in this example 
are: 


K, p = 46900 — kCal - = 3.26 - 1 0 
mole ■ nm' 

K* = 63 — kCCtl •— — 4.38 • 1 O' 10 
mole ■ rad' 


nJ 


bond ■ nm' 

nJ 

angle ■ rad 2 


and the equilibrium bond length, P, , is 0. 1 40 nm. 


(5) 
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4.2 Truss model 


Since it is difficult to express the mechanical strain energy. A', of the truss model in terms of the 
variable truss joint angles that are specified in molecular mechanics (0,-0,), the RVE has been 
instead modeled with extra rods between nearly adjacent joints to represent the interaction 
between the corresponding carbon atoms (Figure 3). In order to represent the chemical model, 
which has bond stretching and variable angles as degrees of freedom, with a truss model that has 
stretching degrees of freedom only, two types of elastic rods, a and b . are incorporated into the 
truss RVE. 

The mechanical strain energy. A 1 , of the discrete truss system in Figure 3 is expressed in the form 
of equation (3): 


Ait \ra ^ -> 


(6) 


where the superscripts correspond to rod types u and b , respectively. Comparing equations (4) 
and (6), it is clear that the bond stretching term in the equation (4) can be related to the first term 
of equation (6) for the rods of type a: 


K* 



(7) 


where it is assumed that p = r‘‘ and P ; = R ] : . However, the second terms in equations (4) and 

(6) cannot be related directly. In order to equate the constants, the chemical bond-angle variation 
must be expressed in terms of the elastic stretching of the truss elements of type b. For 
simplicity, it may be assumed that the prescribed loading conditions consist of small, elastic 
deformations only. This assumption is not an over-simplification for the graphene sheet since the 
deformations for highly stiff linear-elastic materials subjected to many practical loading 
conditions are quite small. 

In order to express the Young's modulus of the rods of type b in terms of the bond-angle 
constant, a simple analysis of the deformation was performed (Figure 4). For small deformations 
of the hexagonal RVE, changes in the bond-angles are small. Therefore, in the second terms of 
equations (4) and (6), it can be assumed that = R‘‘ and that small angle approximations are 
valid. With these assumptions it can be shown that: 


* r>^ 

e.-e.«i~ R 


2R a 


( 8 ) 


Substitution of equation (8) into equations (4) and (6) results in the following approximation: 
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( 9 ) 




Therefore, the Young’s moduli of the two rod types are: 


2 K?R? yh= 3K* 

a i ; ' 2« 


( 10 ) 


The strain energy of the truss model may be expressed as: 


A ’ I 


(ID 


Thus, the strain energy of the truss model is expressed terms of the vibrational potential energy 
constants. 


4.3 Equivalent-plate model 

Working with the assumptions discussed herein, the next step in linking the molecular and 
continuum models is to replace the equivalent-truss model with an equivalent-continuous plate 
with a finite thickness (Figure 3). This replacement is accomplished by determining and 
equating the strain energies of the truss and continuum models for a specific set of applied loads. 
For the case of the RVE shown in Figure 2, a direct solution of strain energy of the equivalent 
plate in terms of displacements is quite difficult to obtain in closed form. A simple and accurate 
way of calculating the strain energy of the truss and continuum RVE representations is to model 
them numerically by using the finite element method [32]. For a given set of loading conditions, 
the strain energies of the two models may then be quickly calculated and compared. Once this is 
accomplished, then the finite thickness of the equivalent plate is determined, and therefore, the 
effective thickness of a graphene sheet is known. 

In order to satisfy' the first condition of the criteria specified in section 3.3, it was assumed that 
the equivalent plate may be modeled using plate-like finite elements. To satisfy the second 
condition, it was assumed that the nodes of the two models are located at the same points, with 
the same degrees of freedom. A simple configuration that satisfies this requirement is shown in 
Figure 3, where the finite element nodes are located at each truss joint on the RVE edges (the 
intersections of the type h rods are not joined) and at the comers of the 4-noded equivalent-plate 
elements. These nodes are allowed to translate in the x, and x 2 directions only (The origin of the 
coordinates is indicated in Figure 3 and located at the centroid of the RVE). Therefore, the 
degrees of freedom of the two models are identical. For condition 3, the displacement gradients 
of the truss elements and the edges of the continuum elements were assumed to be linear. 
Finally, to satisfy condition 4, the two finite element models were subjected to identical loading 
conditions with the total strain energies calculated based on the nodal displacements. This step 
was performed iteratively while optimizing the plate thickness, which was the only available 
adjustable parameter. The corresponding plate thickness is the assumed effective thickness of 
the graphene sheet. 
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While the mechanical properties of the truss elements have been determined as described above, 
those of the graphene sheet were taken from the literature. Reliable values of the in-plane 
mechanical properties of graphite have been known for quite some time since they can be 
measured macroscopically, i.e., without any assumptions regarding the graphene sheet thickness. 
For this example, typical values of the mechanical properties were used [33-35]: 


Y K =1030GPa 
v* =0.170 


where T 8 and are Young's modulus and Poisson's ratio of graphite, respectively. 


4.4 Boundary conditions 

In order to determine an effective plate thickness, the truss and continuum models were 
subjected to three sets of identical loading conditions for each strain energy calculation. For 
each set of loading conditions, a corresponding effective thickness was determined. The loading 
conditions correspond to the three fundamental in-plane deformations of a plate, that is. uniform 
axial tension along x\ and jo and pure shear loading. It was assumed that these loading 
conditions were applied to a macroscopic graphene sheet, i.e. a graphite sheet that has 
dimensions on the macroscopic scale so that there are a very large number of RVEs contained 
within the entire sheet. In order to numerically calculate the strain energy of the RVE under 
these loading conditions, the corresponding periodic boundary conditions for the these 
fundamental loadings were determined and applied to the macroscopic graphene sheet. 


Using the method developed by Whitcomb et al. [36] and Chapman and Whitcomb [37] for 
woven composites, the boundary conditions for the graphene RVE (Figure 2) was determined in 
terms of the overall macroscopic displacement field. The periodic displacement conditions for 
the graphene RVE are given by: 


«r(*.+dr)=»n; 


dr, 


(13) 


where u[ a> is the displacement of a specific RVE face; d[ a> is the vector of periodicity; u k is the 
macroscopic, volume-averaged displacement; x, is the RVE coordinate system located at the 

centroid of the RVE (Figure 3); the superscript a is the RVE face under consideration; and the 
subscripts k,l,m = 1,2, where repeated indices imply summation over the range of the index. The 
vector of periodicity is any vector that connects two equivalent points in a periodic array. The 
vectors of periodicity of particular interest here are those that connect points in adjacent RVEs. 
Figure 5 shows the vectors of periodicity for each set of faces of the RVE. It should be noted 
that dk 1 1 is parallel to a 2 and dk is parallel to r &\ (see Figures 1 and 5). The vectors of periodicity 
for each face are defined with respect to the x* coordinate system by: 
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( 14 ) 


di" = 


i(2) _ 

u k — 


r R‘‘yJ 3 

2 ' 2 


Rfyfe 3 R,‘ 


2 2 

V 

j;-’ =(/?;' V 3 . 0 ) 


Substitution of equation (14) and the geometry of the RVE into equation (13) gives the following 
constraints that represent the periodic boundary' conditions for each group ot opposing faces in 
the graphene RVE: 


.0 


W 3 


-V3x 2 = <" (-y[3x 2 -y[3R%x 2 ) + - 


( 2 ) 


u 


«?' 


\f3x-, + 
’ R;'yJ 3 


2 




3 R“ 

A- 1 . 

r ] 

' 2 

[3x, J 

1 

2 

{ dx 2 J 


-.x, +- 


3r 


= 3x 2 +\J3R‘‘,x 2 j— - 


3 

( du k ) 

+ 3^r 

f afi* \ 

‘ 2 

V ^*1 ) 

2 

\ ^ X 2 J 


(15) 


\ 

x 2 

= I/{ 3) 

3 1 

+ <V3 

M.'i 



l 2 'J 

\ 

J 


Three different sets of displacement boundary conditions were applied to the RVE in this study. 
Each set is discussed below. 


4.4.1 Case 1: uniaxial displacements along ( n/,0 ) 

A macroscopic graphene sheet is subjected to uniaxial loading conditions as shown in Figure 6 . 
An axial displacement A is prescribed on the upper edge while the lower edge is constrained. 
The displacement is relatively small and the h x w plate is perfectly elastic. If it is assumed that 
the global displacements w, and w 2 vary linearly over the global coordinates x, and x 2 , then these 
displacements may be written as: 


(-V ) = V + B 2 

«, ( x* j 5 ,x 2 "t 


(16) 


where B\, B 2 , B 3 , and B 4 are constants which depend on the geometry and applied global 
displacements. The boundary conditions for this case are: 
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(17) 


f 






X P7 

V 2 j 
f 


= A 


x,,-- | = 0 


« 2 

w, ( 0 .x 2 ) = 0 

Since small elastic displacements are assumed, the strains and constitutive equations are [38]: 


e, = 


dE 


8, = 


dll, 

EE 


(18) 


e 'l = ]7f(°.I - V ^) £22 


where £n and 822 are normal strains, and On, and O 22 are normal stresses in the plate. 
Substitution of equations (17) and (18) into equation (16) and solving for the constants B\. ZT , 
Z? 3 . and Z ?4 results in the following displacement functions: 


«■(**) = 

«2 (** ) = 


v"A . 
— x. 


A . A 
— x, + — 
h - 2 


(19) 


The global and RVE coordinates are related by: 

x k = x k + .4 


( 20 ) 


where 4 is a vector indicating the relative position of the two coordinate systems with respect to 
each other. Using equations (15), (19), and (20). the RVE boundary displacements are 
determined to be: 
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V3x 2 +^5.x : + ^1 = m f " * (V3 x 2 +V3/C.* 2 ) + 


V 


: rV3A 


2/7 


W* 2 * 


A /- /?;V3 3/?;' ^ 
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4. 4.2 Case II: uniaxial displacements along (Omi) 

The same macroscopic graphene sheet studied in Case I is now subjected to a different set of 
displacements, as shown in Figure 7. The boundary conditions for this case are: 
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( 22 ) 


If the same assumptions are made for this case as for Case I, then equations (16), (18). and (22) 
may be used to determine the following global displacements: 


.... A „ A 

lV k> w 1 2 

- / . x v'-’A „ 
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(23) 


Using equations (15), (20). and (23). the RYE boundary displacements are determined to be: 
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4.4.3 Case III: (pure -shear displacements) 

The macroscopic graphene sheet is now subjected to pure shear, as indicated in Figure 8. The 
constraints imposed for in-plane shear loading assume that the shearing faces normal to the 
direction of the applied shear are displaced by identical amounts. The displacement is relatively 
small and the h * vr plate is perfectly elastic. If it is assumed that the global displacements ii t and 
i/ 2 vary linearly over x t and x 2 . then these displacements may be written as: 


M, (x k ) = B y x t + B ( x 2 + B 7 
u 2 j -V ) B ^ Aj -f- By \ 2 4“ B w 

The boundary conditions are: 



Substitution of equation (26) into equation (25) results in: 


(25) 


(26) 
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( 27 ) 
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Using equations (15). (20), and (27), the RYE boundary displacements may be determined: 
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To simplify the application of equations (21), (24), and (28) to finite element boundary' 
conditions, without loss of generality , the following assumption can be made: 




(29) 


4.5 Finite element modeling 

The truss and continuum representations of the RVE were modeled using ANSYS* 5.4 [39], In 
the truss model, each pin-jointed extensible rod was modeled using a finite truss element 
(LINK1) with two degrees of freedom at each node (displacements parallel to xi and X 2 ). Rod 
types a and b were assumed to have the same cross-sectional area and different Young's moduli. 
The entire RVE was modeled with the cross-sectional areas of the type a rods divided by a factor 
of 2, since these rods are sharing their total area with adjacent RVEs. The equivalent-plate RVE 
was modeled by using two plane-stress 4-noded quadrilateral elements (PLANE42) with linear 
displacement fields on the edges (Figure 3). The boundary conditions described above were 
applied to each node for a macroscopic graphene sheet height, width, and displacement of lxl 0 (> 
nm, lxlO 6 nm, and 100 nm (corresponding to global uniaxial and shear strains of 0.01%), 
respectively. For both models, the total strain energy of each element was calculated for the 
given boundary conditions, then summed to obtain the strain energy of the entire RVE model. A 
systematic variation of thickness was used to calculate the strain energy of the equivalent plate. 
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and an iterative process was used to find the critical thickness (accurate to 0.01 nm) that resulted 
in the equivalence of the strain energies of the truss and equivalent-plate models. 

5. Discussion 

Using the method described herein, effective thicknesses were determined for the three loading 
cases. For loading cases 1 and II, the resulting effective thickness was calculated to be 0.28 nm. 
For load case III, the effective thickness was 0.24 nm. These values are close to. but smaller 
than, the widely accepted value of the graphitic interatomic spacing, 0.34 nm, and much larger 
than the value suggested by Yakobson el al. [15], 0.066 nm. Even though there is a difference in 
the effective plate thickness based on different loading conditions, the difference is very small 
compared to the range of values currently used. 

It may be assumed that the Young's modulus of a carbon nanotube can be calculated using: 

Y‘"°c(A'"y' (30) 


where A'" is the cross-sectional area of the hollow continuum cylinder with a constant mid-plane 
radius, r . The inner radius, r” , and outer radius, r m , of the tube are: 


» , l 

r = r — 

2 (31) 

m / I 

r = r + — 

2 

where / is the wall thickness. The cross-sectional area of the hollow continuum cylinder is: 

A" 1 = 2m r (32) 

The calculated cross-sectional areas are shown in Table 1. The percentage difference in the 
calculated Young's modulus based on the effective thickness obtained for the three load cases 
with respect to the interatomic spacing is also shown in Table 1. The results suggest that 
measured and calculated values of physical and mechanical properties of carbon nanotubes that 
are dependent on the dimensions of the continuum tube may differ significantly based on the 
assumed geometry. 

The proposed method for continuum modeling may be applied to other nano-structured materials 
such as nanotube-reinforced polymer composites and will yield satisfactory results if the right 
assumptions are made. If the material is subjected to large deformations such that non-linear 
deformation and failure are key issues, then the particular approximations made in this study 
may not be appropriate. 
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6. Conclusions 


A method has been presented for modeling structure-property relationships of nano-structured 
materials. This method serves to link computational chemistry, which is used to predict 
molecular properties, and solid mechanics, which describes macroscopic mechanical behavior 
based on macroscopic material properties. This link is established by replacing discrete 
molecular structures with equivalent-continuum models. It has been shown that this replacement 
may be accomplished by equating the vibrational potential energy of nano-structured materials 
with the mechanical strain energy of representative truss and continuum models. 

As an important example with direct applications to the development and characterization of 
carbon nanotubes, the model has been applied to determine the effective geometry ot a graphene 
sheet. A representative volume element (RVE) of the chemical structure of graphene has been 
substituted with RVEs of equivalent-truss and equivalent-continuum models. As a result, an 
effective thickness of the continuum model has been determined. This effective thickness has 
been shown to be similar to, but slightly smaller than, the interatomic spacing of graphite. 
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continuum 
wall thickness 
[nm] 

continuum 
cross-sectional area 
[nm 2 ] 

% difference of 
Young’s modulus 

[%] 

load cases 1 and II 

0.28 

1.76r' 

21 

load case III 

0.24 

1 .5 1 r* 

42 

average of I, II. and III 

0.27 

1.70r' 

26 

interatomic spacing 

0.34 

2. 1 4r' 

- 


Table 1 - Percent difference in calculated Young's modulus for different wall thicknesses with 
respect to the Young’s modulus calculated with the interatomic spacing. 



Figure 1 - Unit vector notation for a graphene sheet 
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